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ABSTRACT 



Results of measurements from several Arctic field programs and numerical 
models show that clouds affect wind stress during the central Arctic winter by 
changing the longwave cooling of the surface and cloud layers. The longwave 
cooling alters the thermodynamic structure of the lower atmosphere which in turn 
affects the efficiency of momentum transfer to the surface. For typical Arctic 
conditions, wind stress is changed by about 40% one hour after a cloud condition 
change, due to changes in both the surface layer stability and surface layer wind 
speed. The actual wind stress effect due to clouds during this time is a function of 
wind speed, thermal wind, atmospheric boundary layer depth, magnitude of 
radiation change, snow age and, sometimes, snow depth. After several hours, 
surface heat fluxes are no longer important, but the structure of the atmosphere has 
been permanently altered. This affects the wind stress by about 10% to 20% during 
certain situations, but can vary depending on the initial atmospheric structure. 
Measurements of these effects show variations in wind stress associated with clouds. 
Operational and research studies of ice and ocean dynamics will benefit from 
consideration of cloud effects on wind stress. 
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I. INTRODUCTION 



A. BACKGROUND 

This doctoral dissertation is the result of an investigation of the effect of 
clouds on surface wind stress in the Arctic. The two physical phenomena which are 
the focus of this research, clouds and surface wind stress, represent very different 
physical processes, but both play important roles in the Arctic marine atmosphere 
(AMA) and their interactions are poorly understood. 

Clouds affect the AMA in many different ways. The focus of this research is 
on those aspects of Arctic clouds which are most relevant to modeling and physical 
understanding of surface momentum flux or wind stress. Wind stress has special 
importance in the Arctic because it is the primary force driving ice movement 
(Thorndike and Colony, 1982). 

The goal of this research is to answer the following questions, which were 
originally posed in the doctoral dissertation research proposal: 

(1) Are clouds important to wind stress in the Arctic? 

(2) What are the situations when clouds are the most important? 

(3) What is the quantified effect of clouds on wind stress? 

Cloud particles do not directly affect wind stress. Latent heat and radiation 
processes within and outside of clouds cause diabatic heating/cooling of the 
atmosphere and surface. The resulting temperature changes affect surface heat 
fluxes and stability in the lower atmosphere, which are crucial factors in determining 
the geostrophic drag coefficient in the Arctic (Overland, 1985; Overland and 
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Davidson, 1992). Therefore, the answers to the above questions are closely related 
to the thermal structure of the atmosphere and snow/ice surface. 

The very stable atmospheric boundary layer (ABL), which often is present 
during clear weather in the central Arctic, tends to suppress turbulence and 
dynamically de-couple the surface from the rest of the AMA. The surface 
temperature is free to drop until a balance is reached between the upward and 
downward longwave radiation, with a small contribution from conductive heat 
transfer to the surface. If ABL clouds exist, the longwave surface cooling is mostly 
counteracted by radiation from the cloud base. The heat originating from the ice 
surface will be turbulently transported to the top of the clouds where it will be lost 
from cloud top longwave radiation. By changing the location of the major heat loss 
from the surface to the top of the ABL, the clouds have fundamentally altered the 
ABL. 

After a change in cloud conditions, the surface of the snow is subjected to an 
energy imbalance. This will cause the snow surface temperature, hereafter referred 
to as surface temperature, to change and generate turbulent heat fluxes which affect 
both surface layer stability and overall ABL stability. After this initial phase, which 
lasts a few hours, surface layer stability is no longer important, but clouds will 
continue to affect wind stress by changing the overall thermodynamic structure of 
the ABL. 

B. MOMENTUM FLUX AT THE AIR/ICE/SEA INTERFACE 

The transfer of momentum between the atmosphere and the ocean directly 
affects both the atmospheric and oceanic boundary layers (OBL). Momentum 
transfer and the associated shear generate turbulent kinetic energy (TKE) in the 
ABL and the OBL. An increase in TKE can enhance mixing and entrainment, thus 
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thickening the ABL and OBL. The mean currents in the upper ocean are also 
driven primarily by the momentum flux from the atmosphere. 

The formation, destruction and movement of sea ice is controlled by heat and 
momentum fluxes. Campbell et al. (1987) showed that Fram Strait ice movement 
and deformation were primarily controlled by the atmospheric momentum flux 
when surface wind speeds are greater than 6 m/s. Ocean currents, internal stresses 
and surface tilt, which also drive ice motion, are related to wind stress fields. 

Formation of sea ice is a direct result of atmospheric factors such as high 
winds and low temperatures. Melting also can be directly caused by the 
atmosphere, but in the Greenland Sea Marginal Ice Zone (MIZ) most melting 
occurs when the ice is forced over warmer water (McPhee et al ., 1987). In this case, 
the atmosphere causes ice melting indirectly by moving the ice. 

C. ATMOSPHERIC FACTORS WHICH AFFECT MOMENTUM FLUXES 

Momentum flux is controlled by the surface wind speed, surface topography 
and atmospheric stability. The relation between wind speed and momentum flux 
can be parameterized by the equation 

<• = p C a U 2 , (1) 

where r is the surface wind stress, p is air density, U is a surface wind speed and C d 
is a drag coefficient which is a function of surface roughness, atmospheric stability 
and height of U. Considerable progress has been made in recent years in 
developing schemes for estimating the value of C d . 

Although C d as a function of ice conditions is now well known, estimates of 
wind stress in the Arctic are only as accurate as the specification of U. The 
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dependence of stress on the square of U makes wind stress estimates particularly 
sensitive to U. Although surface wind data are typically not available in the Arctic, 
there are enough buoys with pressure measurements to enable the estimation of the 
geostrophic wind speed, U G , and direction with as much accuracy as mid-latitude 

locations (Moritz, 1985). 

The scalar relation between surface wind stress and geostrophic wind can be 
parameterized by use of a geostrophic drag coefficient, C G , which is defined as 



u* 




where u. , the friction velocity, is defined as 



1 




This is similar to C d 1/2 , except that the geostrophic wind rather than the 
surface wind is used. Unlike C d , C G is not easy to determine accurately. Many 
boundary layer effects such as stability, baroclinicity, isobar curvature, depth of the 
boundary layer and horizontal changes of these quantities, as well as surface 
conditions influence C G . Therefore, virtually any phenomenon which affects the 
boundary layer will also affect the value of C G . 



D. CLOUDS AND LONGWAVE RADIATION 

One phenomenon which affects the ABL and C G is clouds. Clouds affect the 

atmosphere in three ways. (1) They release or store latent heat associated with 
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water phase changes; (2) they strongly influence long and short wave radiation 
processes; and (3) they may create hydrometeors (liquid or ice particles large 
enough to be affected by gravity). In the Arctic, liquid and/or ice clouds can exist in 
the boundary layer. This study will analyze only dark season situations with no 
shortwave radiation effects considered. Latent heat effects are minor in the cold 
Arctic winters. Hydrometeors and suspended particulate matter affect radiation 
and surface energy balance, but we do not have good measurements of these 
particles. Therefore, the focus of this study will be on how the longwave radiation 
characteristics of clouds affect wind stress. 

Longwave radiation is affected primarily by (1) water vapor, (2) clouds, (3) 
CO, and other trace gases and (4) aerosol. The amount of water vapor in the Arctic 

atmosphere is limited by the cold temperatures. The other gases are constant, 
except on climatic time scales. Therefore, the effect of changing cloud cover on 
longwave radiation is particularly important in the Arctic. The longwave radiation 
emitted from particles near the top of clouds is greater than the amount absorbed, 
particularly when the air above the cloud is clear. Usually the bottom of a cloud will 
have radiation flux convergence, since the radiation sources below the cloud are 
usually warmer than the cloud base. Both these effects will tend to destabilize the 
cloud layer and lead to larger TKE and entrainment. 

The net effect of the longwave radiation from low clouds is to cool the lower 
atmosphere. Cooling of the lower atmosphere can also occur indirectly due to 
turbulent heat fluxes from a surface cooled by radiation. This is common on clear 
nights over land or non-moist ice. Open ocean and moist ice or snow surfaces have 
mixing (open ocean) or liquid/solid phase changes (moist ice or snow) which 
prevent radiation processes from causing large temperature changes of the surface 
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material. Therefore, in the ABL of the AM A, the major loss of heat due to 
longwave radiation is caused directly by clouds and indirectly by surface cooling. 
The latter is only important when the surface is dry (below freezing) and no clouds 
are present in the ABL. 

This dissertation is arranged as follows. After this introduction, the current 
state of knowledge concerning wind stress and clouds in the Arctic is reviewed. The 
discussion of results begins with an analysis of the effect of clouds on the thermal 
structure of the AMA using experimental data. Next, the short term effects of 
clouds on surface layer stability are analyzed. Following, longer term and complete 
ABL physics are modeled in order to determine the effect of clouds on wind stress 
for realistic situations. Then, statistical relationships between clouds and wind stress 
are examined. The conclusion will summarize the significance of the results and 
potential applications. 
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II. PREVIOUS RESEARCH 



This chapter will examine the state of knowledge concerning wind stress and 
clouds in the Arctic previous to this dissertation. This chapter is divided into three 
sections (1) Results from wind stress studies before the Marginal Ice Zone 
Experiment (MIZEX) programs, (2) Cloud studies based on programs before 
MIZEX and (3) Wind stress and cloud studies from MIZEX and later programs. 
Data from the MIZEX program and later studies were obtained by the investigators 
at the Naval Postgraduate School. 

A. WIND STRESS MEASUREMENTS AND MODELS BEFORE MIZEX 

A summary of results from all wind stress measurements over sea ice reported 
in Western literature before 1983 is available from Overland (1985). Some of these 
results will be mentioned, but the reader is referred to this source for references to 
all studies of C d and summaries of specific values obtained as a result of these 

studies. 

1. Early Studies 

Since the earliest explorations of the Arctic, there has been interest in 
the effect of wind on ice movement (Nansen, 1902; Sverdrup, 1933; Shuleikin, 1938). 
These earlier studies estimated magnitude ratios and angles between wind velocity 
and ice drift. 

The recent emphasis has been on measurement and modeling of the 
ratio between wind speed and momentum transfer. The first estimates of surface 
drag coefficients over sea ice were based on measurements of surface wind speed 
profiles from towers (Untersteiner and Badgley, 1965, Doronin, 1969, Ling and 
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Untersteiner, 1974). Smith et al. (1970) is the first publication reporting C d 
estimates based on direct eddy-correlation measurements using a sonic 
anemometer. 

2. AIDJEX 

The first extensive program to relate wind stress to surface and 
atmospheric factors in the Arctic occurred during the Arctic Ice Dynamics Joint 
Experiment (AIDJEX), which took place in the Beaufort Sea in 1971-1976. Results 
of surface eddy correlation measurements during AIDJEX were reported by Banke 
and Smith (1973) and Banke et al. (1976, 1980), while results from dissipation 
measurements were reported by Leavitt (1980). The value of C d from these studies 

“3 “3 

ranged from 1.14 x 10 to 1.90 x 10 . The locations for these measurements were 
chosen to avoid any wake effects from ice ridges. Therefore, the C d values were too 

low to be used for area averaged wind stress estimates. 

AIDJEX aircraft eddy correlation and momentum integral 
measurements reported by Brown (1977), Katz (1979,1980) and Carsey (1980) 
range from 1.7 x 10* 3 to 2.8 x 10’ 3 . The higher values obtained by the latter methods 
are in good agreement with hindcast estimates of C d from floe trajectories in Coon 

(1980), Hibler (1979) and Neralla et al. (1980). 

3. Other Studies 

In the late 1970's and early 1980’s more surface and aircraft based 
measurements of wind stress over sea ice were obtained using a variety of methods. 
These are reviewed by Overland (1985) and will not be discussed individually here. 
Overland's summary of surface drag coefficients is shown in Table 1. Higher drag 



8 



TABLE 1 

COMPOSITE SURFACE DRAG COEFFICIENTS AS A FUNCTION OF ICE 
AND METEOROLOGICAL REGIME 



Ice Regime 


Characteristics 


7„ -n* 


T„ < -5°. 
Z, < 300 m 


r„ < -5°. 
Z, > 400 m 


Smooth ice 


large, flat floes 


1.5" 


1.5'’ 




Arctic pack 


large range of floe sizes, large 
pressure ridges, C f >0.9 


1.7' 


2.6'' r 




Marginal seas 


broken, first-year ice. C f = 
0.9. occasional big floes 


2.2' 


2.7' 


3.0'’ 


Inner MI7, 


small floes, rafted. Ci - 
0 . 8 - 0 . 9 


2.6' 


3.0' 


3.7 1 


Outer MIZ 


C , = 0.4 

C f = 0.5, rubble field 


2.2' 

2.R’” 







Footnotes refer to the data sets that form the bnsis for the coefficient. Numbers represent a 
subjective median for a range of values. Frror estimates are not possible. 
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coefficients were obtained over sea ice in the marginal seas than in the central 
Arctic. The highest values were found in the MIZs. 

4. Wind Drag Coefficient Models 

Arya (1973,1975) partitioned the surface momentum stress into a skin 
drag and a form drag. The skin drag was caused by the regular floe surface while 
the form drag was due to pressure ridges. The surface measurements from AIDJEX 
were considered to be measurements of skin drag only. Guest and Davidson 
(1984b) extended this model to include the effect of floe edges, which are usually 
the dominant roughness elements in the MIZ. 

Banke et al. (1980) related C d to surface roughness elements after 

filtering wavelengths greater than 13 cm. Brown (1981) modeled the geostrophic 
drag coefficient using a two layer ABL model containing a surface layer and a 
modified Ekman layer. His model included the effect of secondary circulations and 
thermal wind. Above a classic surface layer, the wind profile was calculated by using 
similarity functions which were empirically determined from AIDJEX 
measurements. The value of C G was sensitive to stability at near-neutral conditions. 

The effect of clouds was not considered. 

B. CLOUD STUDIES 

The effect of clouds on the ABL has been the object of a considerable 
research in recent years. Stull (1988) provides a good review of our basic knowledge 
concerning clouds and the ABL. Almost all observational cloud studies have 
occurred in the mid-latitudes and tropics. Many of the observed effects of clouds on 
the dynamics of the ABL can be applied to the Arctic, although differences in 
temperature, cloud condensation nuclei (CCN) size spectra and other factors must 
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be considered. The most common type of clouds in the Arctic are stratus, 
stratocumulus and fog (fog will be considered to be a type of cloud). These types of 
clouds are the focus of this dissertation. 

1. Early Studies 

Vowinckel and Orvig (1970) compiled cloud statistics for the marine 
Arctic which showed a high percentage of stratus in the summer. Jayaweera and 
Ohtake (1973) observed that Arctic stratus usually occurred in several layers, unlike 
mid-latitude stratus. This layering was modeled by Herman and Goody (1976). 
They show that the layering can be caused by shortwave radiation which is trapped 
inside a cloud that is opaque to longwave radiation. They propose that diurnal 
effects prevent similar layering in mid-latitude clouds. 

2. AIDJEX 

The stratus cloud coverage during AIDJEX was anomalously low in 
1975, due to unusual northerly winds (Jayaweera, 1977). Jayaweera also describes 
results from May 1976 Cessna 180 aircraft missions which found that the inversion 
would lower to the center of the stratus clouds after a few days. 

AIDJEX included two radiation missions by the NCAR Electra aircraft 
described by Herman (1977). He estimated Arctic stratus shortwave radiation 
parameters such as bulk values for reflectance, transmittance and absorbance, and 
other parameters such as single scattering albedo, absorption optical depth and 
multiple scattering parameters. 

Using data from the same flights, Herman (1980) determined various 
longwave parameters for the Arctic stratus, including mass absorption coefficients 
for selected liquid water distributions. He concluded that the longwave radiative 
properties of Arctic stratus are similar to mid-latitude stratus. 
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3. Arctic Stratus Cloud Experiment 



During June 1980, six NCAR Electra missions were flown over the 
Beaufort Sea as part of the Arctic Stratus Cloud experiment (ASC). This was the 
most comprehensive program to specifically study Arctic clouds. The first direct 
measurements of cloud liquid water and droplet spectra were made. There were no 
surface based measurements during ASC. 

Tsay and Jayaweera (1984) described the cloud morphology and drop 
size spectra in the clouds. Tsay and Jayaweera (1983) emphasized the different 
radiative properties resulting from wide variations in droplet spectra. Herman and 
Curry (1984) analyzed the effect of the clouds on shortwave radiation with the aid of 
a theoretical model. Curry and Herman (1985a) analyzed the longwave radiation 
properties. These papers contain detailed tables and figures of cloud characteristics 
and radiation measurements. There was a large variability observed in the cloud 
characteristics on different days so that a "typical" summer Arctic stratus cloud could 
not be described. Large-scale factors which affect the occurrence of Arctic stratus 
were examined by Curry and Herman (1985b). 

Curry (1986) examined the interactions between turbulence, radiation 
and microphysics in Arctic stratus by examining four case studies of the ASC flights. 
She has several conclusions which are important to understanding how Arctic clouds 
might affect surface wind stress. (1) The cloud layer is often de-coupled from the 
surface, and several cloud and/or fog layers may exist. Therefore the entire ABL 
cannot be modeled as a mixed layer. (2) The clouds themselves are well mixed due 
to cloud top radiative cooling and latent heat but not due to surface fluxes. (3) The 
longwave cooling in the clouds always exceeds the warming from shortwave 
absorption. The magnitude and vertical distribution of cooling is sensitive to size 
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spectra and liquid water amount at the top of the clouds. (4) The direct effects of 
entrainment on temperature and liquid water are not obvious lower than 50 meters 
below the cloud top. (5) There is no evidence of inhomogeneous mixing as 
described by Baker and Latham (1979). (6) Surface fluxes of moisture do not 

contribute significantly to maintaining the clouds. Gravitational settling has a strong 
effect on the particle size spectra in different regions of the clouds. (7) The spectral 
dispersions of droplet radii are very large due to a variety of mechanisms. (8) Only 
10% of the cloud-top cooling is balanced by turbulent sensible and latent fluxes 
from below. The rest of the heat comes from entrainment and droplet growth. 
Curry does not mention subsidence. 

4. Other Cloud Studies 

A major current program to study marine stratocumulus clouds is called 
the First ISCCP Regional Experiment (FIRE; ISCCP = International Satellite 
Cloud Climatology Project). Much of the research related to this program is based 
on a field program which occurred off the coast of southern California during the 
summer of 1987 (Kloessel et al., 1988). 

Borisenkov et al. (1985) numerically modeled the influence of Arctic 
clouds on large scale meteorology. They found Arctic stratus clouds would increase 
precipitation in Europe and other places far from the Arctic. 

C. WIND STRESS AND CLOUDS: MEASUREMENTS AND MODELS SINCE 

MIZEX 

1. Surface Wind Stress Measurements 

Monin-Obukhov surface layer similarity theory (reviewed by Dyer, 1974) 
applies to surfaces that are horizontally homogeneous. But recent evidence (see 
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Stull, 1988) shows that the theory can also be used when some horizontal variability 
is present. Therefore, more recent Arctic surface measurements were obtained with 
ridges, ice floe edges, or other roughness elements upwind, instead of in the center 
of large flat floes as in earlier studies. 

During the Marginal Ice Zone Experiments of 1983 and 1984, (MIZEX- 
83 and MIZEX-84) many geophysical studies were carried out (MIZEX Group, 
1986). Results from wind stress measurements were reported by Guest and 
Davidson (1984a, b, c, d, 1985, 1987a), Davidson and Guest (1986,1987), Davidson 
and Geernaert (1984a, b, 1985), Fairall and Markson, (1987) and Anderson, (1987). 
All the studies show that C d increased with ice concentration in the MIZ. Guest and 

Davidson (1987a) and Anderson (1987) measured higher drag coefficients over 
rough ice than had been previously reported. The highest MIZ values of C d , 5-6 x 

10' 3 , were measured in regions of very rough ice that had been broken up and rafted 
due to swell action. 

During the MIZEX-87 spring field program (MIZEX '87 Group, 1989) 
the air temperature was well below freezing and new ice was forming. Many stages 
of ice formation were encountered. Guest et al. (1988), Davidson and Guest (1988) 
and Guest and Davidson (1991a) determined values of C d as a function of all ice 

types (Table 2). 

The highest values of C d ever measured over sea ice occurred during the 
Coordinated Eastern Arctic Experiment (CEAREX) drift phase. Multi-year ice had 
undergone extensive deformation as it was pushed, along with the vessel 
Polarbjoem, toward Kvitoya island. Pressure ridges as high as 4 m were common 
around the ship. An average C d value of 8 x 10’ 3 was obtained during a two-day 
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TABLE 2 

THE ROUGHNESS LENGTH, z ( „ AND NEUTRAL DRAG COEFFICIENT, C dn 
FOR VARIOUS ICE AND SEA SURFACES 



Q„ x io 5 



.Ice Type 


x I0 1 m, Median 


Median 


Minimum 


Maximum 


Grease 


0 0027 


0.7 


0.6 


1.1 


Nilas 


0.45 


1.6 


1.4 


1.9 


Pancake 


Diameter <0.75 m 


0.016 


0.9 


0.7 


1.3 


Diameter 0.75-1.5 m 


0.45 


1.6 


1.1 


2.2 


Diameter >1.5 m 


2.8 


2.4 


1.9 


2.9 


Fused 


1.0 


1.9 


1.5 


2.6 


Young 


Smooth 


2.4 


2.3 


1.9 


2.7 


Rough 


7.5 


3.1 


2.6 


3.6 


First year 


Very smooth* 


0.33 


1.5 


1.2 


1.9 


Smooth 


13 


2.0 


1.6 


2.4 


Rough 


7.5 


3.1 


2.2 


4 0 


Very rough 


21.0 


4.2 


3.1 


5.0 


Multiyear 


Very smooth* 


0 33 


1.5 


12 


1.9 


Smooth 


2 0 


2.2 


1.9 


2.5 


Rough 


10 0 


3.4 


2.5 


4.1 


Very rough 


27.0 


4.6 


3.6 


5.5 


Fxtremely rough 


110.0 


8.0 


6.7 


9.1 


Ice free water 
(Steady state open ocean)t 


U = 10 ms ~ 


0 15 


1.3 


na 


na 


IJ = 25 rm ~ ' 


1.3 


2.0 


na 


na 


Bast Greenland Sea 
( V < 12 m < ') 


All wind directions 


0 80 


1.8 


0.7 


3.0 


Ice upwind 2-10 km 


0.23 


1.4 


1.1 


1.8 



Valuer arc based on author's measurements unless otherwise noted and are only from periods when 
lowest inversion was higher than 150 m. At least 85% of the stated ice type was upwind. 

•Includes studies summarized by Overland 1 1985). 
tFrom Smith (1988). No ranges given. 



15 



period. These recent measurements show that C d can have higher values than 
reported by Overland (1985) in certain regions. 

As a result of these studies, the value of the surface drag coefficient can 
be accurately estimated for variations in ice condition. For this to be useful, 
methods must be developed to identify types and concentrations of ice in the Arctic. 
The ice maps currently produced by NOAA give only general descriptions of ice 
types and concentrations. 

2. Observations of Clouds During MIZEX and CEAREX 

During the MIZEX or CEAREX programs continuous surface 
measurements, observations and rawinsonde profiles provide a large data base to 
study many aspects of clouds in the Arctic (Davidson et al, 1984; Lindsay, 1985; 
Guest and Davidson, 1988). During certain periods, there were SODAR 
measurements, aerosol measurements and aircraft missions with cloud physics 
instrumentation. 

Guest (1985) developed techniques for forecasting fog in the MIZ. 
Stratus or low fog is virtually always present during on-ice winds in the MIZ and fog 
forms from stratus lowering. With parallel wind flow, boundary layer fronts are 
common, and the onset of fog or stratus will be abrupt. 

Statistical studies by Guest and Davidson (1987b, 1988) show that the 
presence and thickness of clouds were correlated with the height of the inversion 
base, Zj. When stratus clouds were present, Z s was higher and was less correlated 

with wind speed than clear sky cases (Figure 1). 

Guest et al (1988) compared median Zj values for on-ice and off-ice 
winds in the MIZ (Figure 2). At all relative locations, Zj medians were higher for 
on-ice winds than off-ice winds. The difference in Zj becomes greatest over the 
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Figure 1 Composite Median Inversion Base Heights, Z-, from MIZEX-84 Divided 
into Wind Speed and Cloud Thickness Categories. Note that during clear (no 
cloud) conditions the inversion was lower and more greatly affected by wind speed 
than during cloudy conditions. 
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Figure 2 Composite Median Z ; Values Across the MIZ for Different Wind 
Directions during M1ZEX-87. Triangles, circles and squares represent off-ice (left- 
to-right), parallel-right (out-of-paper) and on-ice wind conditions, respectively. 
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pack ice. During on-ice winds, fog or stratus inevitably would exist over the ice. 
The destabilizing effect of the clouds may be a major reason for the observed 
difference in Z ; between wind regimes. (There are also influences due to advection 
and subsidence.) If the clouds affect Z i( they will also affect surface momentum flux, 
because a higher Zj is associated with more efficient transport of momentum to the 
surface. 

3. Boundary Layer Models 

Although surface drag coefficients are well known and straightforward to 
use, their use requires measurements of surface wind speed. Estimating wind stress 
from a pressure field is much more complicated. This involves the specification of 
the geostrophic drag coefficient, C G , which requires some type of ABL model. 

Overland (1985) used a steady-state one-dimensional ABL model to 
determine C G and the turning angle between the geostrophic and surface wind, a, as 

a function of a mechanical mixing stability parameter. He found that when the 
atmosphere became very stable, the surface wind tends to de-couple from the upper 
winds, and C G is no longer affected by surface roughness. These results are 

applicable in the Arctic pack ice away from the MIZs. Clouds effects were not 
considered. 

In the 1980's, the focus of Arctic studies was on MIZs. The few ABL 
models for the MIZ that have been published will be reviewed. Overland et al. 
(1983) used a primitive equation slab model to explain observations of the Bering 
sea MIZ in March during off-ice winds. Reynolds (1984) examined the same case 
with a similar model, but included radiative cooling from clouds. He attributed the 
observed variations in temperature, wind and ABL height to drag variation over the 
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ice and heat flux over the open ocean. Clouds and heat flux through the ice were 
not important. 

Andreas et al. (1984) described a MIZ ABL based on a rawinsonde 
cross-section of the Antarctic MIZ during off-ice winds in the Weddell sea. Bennett 
and Hunkins (1986) simulated this case with a two-dimensional, multi-level model 
which includes radiation. They concluded that adiabatic lifting due to convergence 
and longwave cooling from clouds dominate the cooling in the ABL. Although 
there was some discussion concerning model details (Andreas, 1987; Bennett and 
Hunkins, 1987), all authors agreed that there were large changes in the wind stress 
across the MIZ due to roughness and cloud effects. 

Chu (1986a, b, c, 1988a, b, c) has modeled several mechanisms by which 
air-ice-sea interactions may cause ice features in the MIZ. Chu et al. (1990) and 
Chu and Garwood (1990, 1991) investigated feedback mechanisms between clouds 
and surface fluxes over ice-free oceans. 

Brown (1986) applied his one-dimensional, two-level ABL model to an 
off-ice wind case in the Fram Strait MIZ. A realistic surface roughness and 
temperature field was used. The greater roughness in the outer MIZ slows and 
backs the surface wind. A greater upward heat flux at the surface has the opposite 
effect. 

Recent results of two dimensional multi-level MIZ ABL models by 
Glendening (1992) and Kantha and Mellor (1989) show the same general results for 
off-ice winds. Kantha and Mellor also examined the MIZ ABL during on-ice, 
parallel-left and parallel-right wind regimes. The above studies show the surface 
stress field is very complicated in the MIZ because of the changes in the ABL due to 
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surface roughness and horizontal temperature variations. Cloud effects have not 
been modeled by these researchers. 

Glendening's model of the MIZ ABL during off-ice winds shows that Z ( 
is controlled by surface wind stress over the ice regions and surface temperatures 
(heat fluxes) over open ocean regions. He has not examined very stable cases when 
surface conditions become de-coupled from the upper ABL. 

Modeling stable ABL situations is more difficult because of the 
intermittent nature of the mixing events. Overland (1988) approaches the problem 
with a one-dimensional, multi-level, TKE-mixing length model which simulates the 
ABL over ice in the winter. When cloud radiation effects are introduced, the 
boundary layer becomes unstable and surface stress increases by approximately 
20 %. 

4. Wind Stress in the MIZ 

The above models, as well as observations (Davidson and Guest, 1988) 
show that the surface wind stress field does not match the surface drag coefficient 
field because there are changes in the surface wind speed across the MIZ. The wind 
speed at the surface is affected by changes in stability and surface roughness. By 
using observations and Overland's (1985), Brown and Liu's (1982) and Brown's 
(1986) ABL models, Campbell et al. (1987) and Guest (1988) showed that the 
combined effects of horizontally varying surface wind speeds and roughness result in 
a wind stress field which is complicated and dependent on the wind direction 
relative to the ice. Typically, the wind stress is greater over the open ocean than the 
pack ice, even though C d has an opposite distribution. This fact has been neglected 

by several MIZ ice movement models (Roed, 1983; Roed and O'Brien, 1983; 
Hakkinen, 1986a, b; Smith et al. 1988). These models use an unrealistic constant 



21 



surface wind speed and direction across the MIZ, which results in greater stress over 
the ice regions. 

Guest (1988) also shows that clouds can have a large effect on surface 
wind stress if they change the stability of the surface layer from stable to unstable. 
This effect of clouds will be further examined. 
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III. AN EXAMINATION OF THE RELATION BETWEEN CLOUDS AND 
ATMOSPHERIC BOUNDARY LAYER THERMODYNAMIC STRUCTURE 



Clouds influence wind stress by their thermodynamic effects on the 
atmosphere. This chapter will examine the thermodynamic effects of clouds, which 
form the basis for understanding the effect of clouds on wind stress. 

Divergence of longwave radiation at cloud tops and/or at the snow surface has 
a significant effect on the thermodynamic structure of the lower atmosphere over 
pack ice during the Arctic winter. Overland and Guest (1991), hereafter OG, 
showed how the temperature at the surface was primarily controlled by longwave 
radiation. Heat conduction through the ice and leads are less important but do 
prevent extremely cold temperatures ( < -40 C) from existing over sea ice. The layer 
of air just above the inversion, or isothermal layer, contains the warmest 
temperatures in the atmosphere (Figure 3). It will tend to experience more 
radiational cooling than the rest of the atmosphere. Warm horizontal advection 
must counteract the longwave cooling in the isothermal layer, on the average. 

The actual location of the maximum longwave cooling in the lower 
atmosphere depends on the cloud conditions. The location of the cooling affects the 
stability and therefore the dynamics of the lower AMA. When conditions are clear, 
the maximum longwave cooling occurs in a thin layer at the top of the snow/ice. 
When clouds exist, the maximum radiational cooling occurs at the top of the upper 
cloud layer. There is a strong correlation between amount of cloud cover and 
downward longwave radiation (Figure 4). When low overcast is present, all of our 
Arctic data indicate a clear trend: the downward longwave radiation at the surface is 
within ± 10 Wm' 2 of the blackbody radiation at cloud bottom. Therefore, it can be 
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Temperature (deg C) 



Figure 3 A Typical Temperature Sounding From the CEAREX Drift. The mixed 
layer, lower 200 m, inversion layer, 200 to 500 m, and the isothermal layer (500 to 
1200 m) are shown. 
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DOUNIJARD LONGUAUE RADIATION <U/M~2> 



Figure 4 Frequency Distribution of Downward Longwave Radiation during the 
CEAREX Drift for Clear (small squares) and Low Overcast (large diamonds) 
Conditions. Note that the peak frequency for cloudy conditions is about 75 Wm 2 
higher than for clear conditions. 
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assumed that all but the thin (less than 100 m) low level cloud layers that exist in the 
Arctic are essentially black to longwave radiation. The examination will be based 
on data obtained during the MIZEX and CEAREX programs. 

A. CEAREX DRIFT MEASUREMENTS 

Most of the observational data used for this study were from measurements 
from the Polarbjoem during the CEAREX drift phase, 15 September to 12 
December, 1988. During the drift phase, the Polarbjoem was moored to a large ice 
floe which drifted in the region North and East of Svalbard. Instrumentation was 
located both on a large ice floe (alpha floe) to which the Polarbjoem was moored, 
and on board. 

A meteorological station was located on a platform extending forward of the 
Polarbjoem 's bow mast at a height of 14 meters above sea level. This station 
measured wind speed and direction, temperature, and relative humidity. 
Temperature data were accurate to within 1.0 C. The humidity measurement was 
not reliable. Wind speeds were accurate to within 0.3 m/s and direction to within 10 
degrees, although during periods of unfavorable wind direction, errors may have 
been larger. 

The radiation system measured both downward infrared and downward solar 
irradiance, and was located on the ice until the floe was crushed on November 15. 
At this time the system was relocated aboard the Polarbjoem. The system remained 
operational until December 12. The sensors would often frost quickly; therefore 
hourly cleanings were performed. Radio transmissions would interfere with the 
measurements. The data were manually edited to eliminate poor measurement 
periods. The instruments were modified so that errors due to temperature 



26 



differences within the sensor could be identified and corrected. The resulting 

2 

irradiances were generally accurate to within 5 watts/m . 

Rawinsondes provided vertical profiles of temperature, humidity, and wind 
speed and direction twice daily throughout the experiment, with additional 
soundings during periods of extreme or unusual weather. The sondes were 
equipped with thermistors that measure temperatures to within 0.2 C. The humidity 
sensor consisted of a specially coated glass plate, the resistance across which varies 
with humidity. It was accurate to about 5%. Sonde wind directions were generally 
accurate to within 20 degrees. 

A 6 meter profile mast, located on an adjacent floe, measured temperature 
and wind speed at 4 levels. This was used to provide wind stress and heat flux 
measuirements when the wind direction was favorable. The heat flux was estimated 
to be accurate to 5 Wnr 2 and the wind stress to 20%. Surface temperature was 
measured with a thermistor placed on the snow surface and was accurate to 1.0 C.. 

B. THERMODYNAMIC COUPLING OF THE SURFACE, CLOUD LAYER 

AND ISOTHERMAL LAYER 

This section will use the results of measurements from the CEAREX drift to 
examine the relationship between surface temperature, air temperature, ABL 
temperature and isothermal layer temperature over ice in the Arctic winter. (In a 
strict sense, the term "air temperature" refers to the temperature at 10 m above the 
surface. In practice, the temperature at the top of the ice profile mast, 5.2 m, and at 
the ship mast, 14 m, were used to determine air temperature.) The differences 
between these temperatures determine the stratification of the lower atmosphere 
which in turn affects momentum transfer to the surface. The effect of clouds on the 
stratification will be examined. 
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The surface temperature and the air temperature are closely linked (Figure 5) 
because surface sensible heat flux quickly counteracts any temperature difference 
between the snow and near-surface air. The one period when there was a difference 
of greater than 2 C, 2-6 November 1988, occurred during very low or zero winds, 
when turbulent surface fluxes were suppressed. Usually, enough turbulence exists to 
support surface layer heat fluxes and therefore the difference between surface 
temperature and ABL potential temperature is usually very small. 

OG showed how surface temperature, and therefore air temperature, is 
thermally coupled to a radiational boundary layer (RBL). The RBL is characterized 
by an isothermal temperature layer above an inversion layer. The isothermal layer 
extends to approximately 1.5 km elevation. The surface is prevented from becoming 
extremely cold because the heat that is lost by upward longwave radiation is 
replaced by downward radiation from the isothermal layer, plus a small amount of 
heat conduction upward through the ice/snow. 

OG did not consider the effect of clouds, which, in the isothermal layer or 
ABL, would greatly increase the downward radiation at the surface compared to 
clear sky conditions. This causes a very close thermodynamic coupling between the 
clouds and the surface. This coupling is illustrated in Figure 6, a time series plot of 
surface air temperature vs. sky temperature. Sky temperature, T sky , is defined as the 

temperature of a blackbody corresponding to the measured downward radiation at 
the surface, R downsfc 
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Figure 5 Time Series of Surface Temperature, 5.2 Meter Air Temperature, Wind 
Speed and Air-Snow Temperature Difference during the CEAREX Drift. 
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Figure 6 Time Series of 14 Meier Air Temperature (dashed line) and Sky 
Temperature (solid line) during the CEAREX Drift. 



30 



where a = 5.67 x 10' 8 Wm' 2 K 4 is the Stefan-Boltzmann constant. The air 
temperature, T air , in Figure 6, is usually very close to the sky temperature, except for 

a few periods of colder sky temperatures. 

The periods when T air was close to T sliy were cloudy periods, while the periods 
when T s ky was colder occurred during clear periods. This is illustrated by 
scatterplots of T air . vs. T sky for overcast and clear sky conditions (Figures 7 and 8). A 
link between T air and ^sky exists for all sky conditions. This is because the dry snow 
surface has a low heat capacity and low thermal conductivity and therefore the snow 
surface temperature quickly responds to changes in surface energy balance caused 
by changes in downward radiation. The snow surface is then closely linked to the 
ABL air temperature by turbulent heat flux as discussed above. 

During low level overcast periods, a better correlation between air 
temperature and surface temperature exists (Figure 8). T air and T s ky are never more 

than a few degrees different from each other. During these overcast periods, the 
surface is not only linked by radiation with the surface, it is also turbulently linked if 
the cloud bottom is within the ABL. ABL turbulence effectively transfers heat 
between the surface and the top of the ABL during overcast conditions. 

A conclusion of OG was that longwave radiational cooling of the snow surface 
causes the lower atmosphere to cool during the dark seasons. The temperature of 
the air near the surface is determined by a balance between upward radiation at the 
snow surface and downward radiation from the isothermal layer (with a small effect 
due to heating from the ocean through leads and ice). This study proposes that 
these conclusions be modified for situations when the ABL is cloud-capped. During 
these periods, important radiational exchanges are occurring at the cloud top rather 
than at the surface. 
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Figure 7 Scatterplot of Sky Temperature vs. 14 Meter Air Temperature for Clear 
Sky Conditions during the CEAREX Drift. The line represents surface temperature 
equal to sky temperature. 
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Figure 8 Same as Figure 7 but for Low Overcast Conditions. 



Overland and Davidson (1992), hereafter referred to as OD, found that wind 
stress was affected by a measure of the external or background atmospheric 
stratification based on the difference in temperature between the surface and an 
upper level (900 or 850 mb), assumed to be in the isothermal layer. When ABL 
clouds exist, the surface is no longer directly connected by radiation to the upper 
level. In this case, an external stability parameter defined by the difference in 
temperature between the cloud top (instead of the surface) and an upper level may 
be appropriate, since it is at the cloud top that the direct radiational linking between 
the surface and the upper levels is actually occurring. 

C. CLOUDS AND BOUNDARY LAYER THERMAL STRUCTURE 

To determine the effect of clouds on wind stress, it is important to understand 
how clouds affect or are associated with ABL structure. Two parameters for 
describing ABL structure were derived from hundreds of different MIZEX and 
CEAREX rawinsonde soundings. These were (1) height of the lowest inversion 
base, Zj, and (2) "type" of boundary layer, described below. 

1. Inversion Height - Average Values 

Low cloud cases have a higher average Z, than clear sky cases. This is 

shown in Table 3 for locations within the pack ice, in the MIZ and in the adjacent 
open ocean. Average inversion heights are considerably higher for cloudy skies 
compared to clear skies. This is true for the various locations shown in Table 3 and 
also is true for any data set grouped by season, platform, wind direction regime (off- 
ice, on-ice) or wind speed regime (not shown). 
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TABLE 3 

MEAN AND MEDIAN Z, FOR DIFFERENT LOW CLOUD CONDITIONS 

AND LOCATIONS 



Cloud Condition 


Pack ice* 


Location 

MIZ 2 


3 

Open Ocean 


Overcast 


median 


220 


400 


610 




mean 


270 


480 


730 




sd 4 


210 


390 


470 


Clear 


median 


100 


60 


420 




mean 


160 


220 


450 




sd 4 


240 


370 


340 



*50 km or more from ice edge (primarily CEAREX drift data) 

? 

<50 km from ice edge but still over sea ice 

3 

over open ocean but within 200 km of ice edge 

4 

sd = standard deviation These values represent natural variations, not 
experimental errors. 



35 



The large standard deviations of the Zj values show that other factors 

besides clouds affect ABL structure. An analysis of variance shows that 22% of the 
variation in Zj values not associated with distance from the ice edge can be 
attributed to cloud conditions. The higher average Zj values associated with low 
clouds are not due entirely to the direct effect of clouds on the ABL. Cloud 
occurrence is associated with warm advection, low pressure systems, high wind 
speeds and moist marine air masses. All of these would be associated with higher Zj 

values even without cloud effects. 

In this section, Z ( vs cloud relationships were examined, and in the 
previous section, the coupling between longwave radiation and ABL air 
temperature for different cloud conditions was discussed. Results from both topics 
can be combined by plotting the difference between T s ky and T air vs. Z ; for clear and 

low overcast conditions (Figure 9). Several differences between overcast and clear 
conditions are apparent. Overcast conditions are characterized by Tsky being nearly 
equal or slightly greater than T air , and by having higher average inversion bases 
when compared with clear cases. Surface-based inversions (Z ; = 0) are extremely 
rare during cloudy skies, but clear skies do not guarantee that a surface inversion 
will be present. 

2. Inversion Height and u. Scaling 

It has been observed since early explorations (Sverdrup, 1933) that the 
value of Zj is related to surface wind speed. Increasing wind speed generates more 

turbulent kinetic energy, TKE, due to turbulent shear production, which is able to 
extend mixing to higher levels. Attempts to scale or predict Z ( or the ABL depth, h, 

for stable or neutral ABLs invariably use parameters related to the mechanical 
production of turbulence: the friction velocity, u., an eddy viscosity coefficient, K, 
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Figure 9 Scatterplot of the Difference Between the 14 Meter Air Temperature and 
the Sky Temperature vs. Inversion Base Height for Clear (squares) and Low 
Overcast (plusses) Sky Conditions. 
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the bulk Richardson number, Ri B , or a wind speed or wind shear parameter. The 
following scales for equilibrium ABL height, h e , have been proposed for stable and 
neutral ABLs. 

Neutral ABL and Free Atmosphere: 



he 




7 



Ekman depth 



h e 




e.g. Panofsky and Dutton (1984) 



Stable Surface Layer, Neutral Free Atmosphere: 



0 v u* 3 






Obukhov length scale 


kg w' e v ' 


/u*L\j 




( f ) 


Zilitinkevich (1972, 1974) 


Ri B U h 2 9 
gK - O 


Hanna (1969), Wetzel (1982) 



Neutral Surface Layer, Stable Free Atmosphere.- 



* , g de v 

h e a ; N“ = Q ~ OG , Kitaigorodskii (1988), 

Kitaigorodskii and Joffre (1988) 



u. 

The parameter ^ was shown by OD to be relevant for scaling the Arctic 



winter boundary layer, when surface fluxes are usually small and an external or 
background stability limits the height of the ABL. OD used data from the 
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CEAREX drift to verify this relationship. A re-analysis of these data confirms that 
Zj was correlated (R 2 = 0.34) with u. or wind speed during the CEAREX drift data 

(Figure 10). A similar correlation (R 2 = 0.21) occurs for the data set containing all 
MIZEX and CEAREX soundings made by our group in the Arctic (Figure 11) when 
all sky condition cases are included. However, when either data set is divided into 
categories depending on cloud conditions, it is found that the wind speed (or u.) 

dependence completely disappears (R 2 < 0.05) for the low overcast cases (Figures 12 
and 13). 

The wind dependence noted by OD and others apparently occurs only 
during clear sky periods. This lack of wind speed dependence for Zj during overcast 

conditions was observed in the MIZ during the summer in MIZEX-84 (Guest and 
Davidson, 1987b) and in the spring during MIZEX-87 (Guest et al . , 1988) as well as 
over pack ice during the fall/winter CEAREX drift (Guest and Davidson, 1991a). 
Therefore, low clouds appear to reduce the relative influences of surface layer 
mechanical production of turbulence on the Z| values. 

The latter conclusion seems to contradict some basic assumptions about 
the factors that control ABL depth in a stable atmosphere. The lack of a wind 
speed vs. Zj correlation could have several explanations. (1) Z i5 as defined here, is 

not always a good measure of ABL height, h. For the modeling studies, the ABL 
height, h, is arbitrarily defined as the height at which TKE becomes 10% of the 
surface value. Clouds may create a sharp inversion which defines the Z| value when 

mixing is actually limited to a lower level which cannot be identified from the 
rawinsonde data. (2) The boundary layer under clouds in the Arctic is usually 
unstable and therefore buoyant production of TKE may dominate over mechanical 
production at times. (3) Cloud top cooling and wind shear, both of which often 
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Figure 10 A plot of Inversion Height vs. Wind Speed during the CEAREX Drift. 
The solid line is the linear regression fit, the inner dashed lines are the 95% 
confidence intervals of the regression and the outer dashed lines are the 95% data 
range limits. 
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Figure 11 Same as Figure 10 blit using All Data from MIZEX-84, MIZEX-87 and 
CEAREX. 



Inusrsion Baser Haight <m) 



CEAREX Drift Low Ouercest Only 




Figure 12 Same as Figure 10 but for Low Overcast Cases Only. 
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Figure 13 Same as Figure 1 1 bul for Low Overcast Cases Only. 
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occur at the top of the ABL, may be more important than surface processes in 
controlling Z ( during overcast conditions. (4) Advection of moisture at certain levels 
may fix Zj at that level regardless of surface conditions. (5) Cloud top entrainment 
instability, CTEI, (Lilly, 1968) is another process not directly related to surface wind 
stress. A necessary condition for CTEI is that the equivalent potential temperature 
of the air in ABL must be greater than the air just above ABL. This condition is 
rare in the Arctic; therefore CTEI cannot be common. 

Although this sub-section concerned the effect of wind stress on an ABL 
structure parameter, Zj, the focus of this study remains the more subtle reverse 

process: the effect of cloud-influenced ABL structure on wind stress. These results 
have illustrated how the surface wind stress influence on boundary layer structure is 
strongly dependent on cloud conditions. Scaling inversion height with u. is not valid 

during low overcast conditions, which provides further evidence that clouds have a 
strong effect on ABL structure. 

3. Types of Atmospheric Boundary Layers 

Another way to describe ABL thermal structure is to classify each 
measured ABL according to a scheme summarized in Table 4, where the number 
assigned each type is arbitrary. For types 0, 3, and 5, the specification of Zj is 

straightforward since these represent "classic" ABLs where a clearly defined 
inversion exists at the surface or above a surface-based well-mixed layer. For the 
other types the Z ; specification is more ambiguous because either no inversion exists 
or weakly stable layers exist below a more obvious inversion. In these cases, Z ; is 
still defined as the lowest inversion base, but probably is not a good measure of 
ABL depth, h. 
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TABLE 4 

TYPES OF ATMOSPHERIC BOUNDARY LAYERS 



Type 

-1 No inversion in lower 2500 m 

0 Strong surface-based inversion (height of top of inversion ,Z top > 200 
m or potential temperature increase within the inversion, (9 jump > 3 C) 

1 Weak sfc-based inversion below mixed layer (Z top < 200 m and Z top > 
100 m and 0 jump < 3 C) 

2 Elevated inversion above stable layer; stable layer not inverted. 

3 Elevated inversion above mixed layer 

4 Miscellaneous category (complicated structure) 

5 Two inversions with mixed layer between 

0.5 Combination of 0 and 5; A strong surface-based inversion below a 

well-mixed layer capped by an upper level inversion. 
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All rawinsonde profiles we have collected in the Arctic during MIZEX 
and CEAREX (over 1300) have been classified according to ABL type and amount 
of low clouds. The percentage of total of profiles for each type is listed in Table 5. 

As indicated previously (Figure 9), surface-based inversions are rare 
during low overcast conditions and common during clear skies. All the low level 
inversions occurred over ice and during clear skies, with the exception of four cloudy 
cases when strong warm advection over cold ice in the MIZ caused a surface-based 
inversion to form. During the CEAREX drift, there were no surface-based 
inversions, Type 0, with overcast skies. This is a significant observation with respect 
to the wind stress because during these very stable, low inversion periods 
atmospheric stability determines how wind stress is related to geostrophic forcing 
(Overland, 1985; and later in this dissertation). During less stable, higher inversion 
periods, changes in surface roughness are more important than atmospheric stability 
in affecting wind stress. Clouds play an important role in destroying or preventing 
the formation of strong surface-based inversions. 

Multiple upper-level inversions (Type 5) were twice as common for 
overcast vs. clear skies. Most of the Type 5 clear cases had high relative humidities 
at the bottom of the upper inversions, which may indicate evidence of previous 
clouds. Arctic clouds can form mixed layers capped by inversions which are de- 
coupled from the ABL. 

When surface-based inversions existed below an upper-level inversion 
(Type 0.5, 12 cases), only one case had clouds associated with the upper inversion. 
Clouds prevent surface inversions from forming even if they are de-coupled from 
the surface. 
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TABLE 5 



A CROSS TABULATION OF BOUNDARY LAYER TYPE BY LOW-LEVEL 
CLOUD CONDITIONS EXPRESSED AS PERCENTAGE OF ALL CEAREX 

DRIFT SOUNDINGS 1 

/ 

Low Cloud Amount 

2 

ABL Type Clear Partial Overcast 



Surfaced-based inversion 
Type 0 or 0.5 


13.3 


4.4 


0.7 


Inversion above mixed layer 
Type 3 or Type 5 


34.8 


11.1 


26.7 


Two inversions 
Type 5 or 0.5 


23.7 


12.6 


17.0 


Unclear inversion location 
Type 1, 2 or 4 


4.4 


2.2 


2.3 



^Based on 135 rawinsonde profiles. Total of all categories is 153%. This is because 
14% of the soundings were Type 0.5, which are both surface-based and two 
inversions and 39% of the soundings were Type 5 which is both mixed surface layer 
and two inversions. 

^Type numbers are defined in Table 4. 
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D. CLOUDS AND SURFACE FLUXES 



One way clouds can affect surface wind stress is by their influence on the 
surface layer stratification. This section is an examination of how observed sky 
conditions and sky condition changes are related to measured surface sensible heat 
fluxes. Turbulent heat flux at the surface determines surface layer stability which 
affects the value of the surface drag coefficient by the well-verified Monin-Obhukov 
similarity theory (Dyer, 1974). Also, surface heat fluxes can affect wind stress by 
influencing entrainment into the ABL. Entrainment brings momentum 
characteristic of upper-level air into the ABL, causing a short term effect on surface 
wind stress. Entrainment also changes the ABL depth, h, which can have a long 
term effect on wind stress. Results from numerical models of these processes are 
presented in Chapter V. 

Thirty-one periods from the CEAREX drift were identified when the wind 
direction was favorable for surface heat flux measurements, meteorological 
conditions remained constant and several consecutive ten minute heat flux values 
were similar. Low cloud cover was correlated with surface sensible heat flux (Figure 

-9 “2 

14). The average heat flux was -5 Wm ' and 1 Wm during clear skies and low 
overcast skies respectively. A better correlation occurred between measured 
sensible heat flux and downward longwave radiation at the surface (Figure 15) or 
net surface radiation (Figure 16), although scatter remains. The zero-flux calm wind 
cases are not shown. All the measured heat fluxes from the CEAREX drift were 
small compared to typical radiative terms and MIZ or lead turbulent heat fluxes. 
However, the largest changes in the effect of stability on wind stress occurred in the 
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Claud Couer (Octaues) 



Figure 14 Surface Sensible Meat Flux (positive upward) as a Function of Cloud 
Cover Octaves from Selected CEAREX Drift Periods. Lines represents statistics 
described in Figure 9 caption. 
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Sensible Heat Flux (Um A 2) 




Figure 15 Surface Sensible Heat Flux as a Function of Downward Longwave 
Radiation (Wm’ 2 ). Very low winds cases are not shown. 
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Sansibl® H«at Flux (Um~2) 




Figure 16 Surface Sensible Heat Flux as a Function of Net Longwave Radiation 
(Wm' 2 ). 
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near-neutral regime (e.g. Brown and Liu, 1982), so even small differences in heat 
flux and surface stability may significantly affect wind stress. 

Most periods represented by these data had quite constant meteorological 
conditions, including cloud cover. Larger sensible heat fluxes were measured 
immediately after changes in cloud conditions. This is best illustrated in a case 
study (Figure 17) of a period from the CEAREX drift. Note the close correlation 
between cloud conditions and surface sensible heat flux for this case. Cloud clearing 
is closely linked to a change from unstable to stable surface conditions while the 
reverse occurs for cloud forming. Several other similar examples exist from the 
CEAREX drift period. 

Sixteen sharp cloud change events (eight totally clear to low overcast changes 
within an hour, eight low overcast to clear changes) were identified during periods 
when surface flux measurements were available. The surface layer of the 
atmosphere changed from unstable to stable in all the clearing cases. The surface 
layer changed from stable to unstable in five of the "clouding" cases. The other 
three clouding cases showed no significant change in heat fluxes or stability. It is 
obvious that surface fluxes are closely related to changes in cloud conditions, 
particularly immediately after clearing events. 

E. FACTORS INFLUENCING SURFACE HEAT FLUX 

Instantaneous and integrated turbulent surface heat fluxes affect the ABL 
structure and, therefore, the ability of the atmosphere to transfer momentum to the 
surface. It has been shown that changes in cloud conditions are often associated 
with changes in the surface turbulent heat flux. But other factors can also induce 
heat flux changes. In order to evaluate the significance of the effect of clouds on 
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wind stress in the central Arctic winter, the effect of downward longwave radiation, 
i.e. cloud conditions, on surface heat fluxes must be compared with other processes 
that can generate a surface heat flux. This section describes a method for classifying 
a non-neutral ABL according to the mechanism driving the surface turbulent heat 
flux. The classification scheme will then be applied to the CEAREX drift data to 
evaluate the relative importance of changes in surface radiation conditions in 
generating surface heat fluxes relative to other mechanisms. 

1. General Classification Scheme 

With no outside forcing or advection, the turbulent surface heat flux 
drives the T air - T sfc difference toward zero over time. As shown in section A, these 

two parameters are closely linked unless no turbulence exists in the surface layer. 

Surface heat flux counteracts some "outside" forcing mechanism creating a 

difference between T air and T s fc- There are two types of outside forcing 

mechanisms: Type R (radiational) ones that act directly to change T sfc and Type A 

(advective) ones that act directly to change T air . In Type R ABLs, the surface is 

controlling the temperature in the ABL through surface fluxes; changes in Tsfc 

precede changes in T ajr . In Type A ABLs, the temperature in the ABL is controlling 

T s fc through surface fluxes; changes in T air precede changes in T s fc- These two 

mechanisms can be distinguished by comparing the signs of the surface heat flux and 
5 T air 

. For Type R they are the same sign while for Type A they are opposite signs. 
This concept is summarized in Table 6. 
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TABLE 6 

TYPES OF TURBULENT HEAT FLUX EVENTS 



Type R - Radiational 

Downward (upward) turbulent surface heat flux 
Cooling (warming) ABL and surface 

Type A - Advective 

Downward (upward) turbulent surface heat flux 
Warming (cooling) ABL and surface 
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Type R is identified as "radiational" because changes in T sfc not caused 

by sensible heat fluxes are closely linked to changes in radiation conditions. The 
only other factor affecting T sfc is conductive heat flux through the snow/ice. This 

can be considered a passive process because snow and ice conditions do not change 
much on the time scales of a few hours being considered here. Type R ABLs are 
common over land where diurnal variations in temperature are large, particularly in 
the summer when solar radiation is intense and the soil is dry. Because advection of 
soil or snow is insignificant, pure Type R ABLs are by definition caused by a one- 
dimensional process. 

Changes in radiation conditions in the Arctic winter are strongly related 
to cloud conditions. These radiational changes are an order of magnitude less than 
diurnal solar radiation changes at lower latitudes, where Type R ABLs are common. 
However, it will be shown that the winter Arctic ABL often has Type R 
characteristics. This is because the temperatures of the dry snow surface and 
shallow ABL quickly adjust to radiation changes. 

Type A ABLs are caused by any process, other than turbulent surface 
heat flux, that changes the temperature within the ABL. Note that this is termed 
"advective" in Table 6. This is because advection is assumed to be the dominant 
non-radiative process controlling the structure of the ABL. Type A ABLs occur 
over open water regions with SST gradients and anywhere following an atmospheric 
frontal passage. 

2. Scale Analysis 

To evaluate the factors which change ABL temperature, consider the 
following equation for the time rate of change of the vertically-integrated potential 
temperature of the ABL. 
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The overbar represents a vertical integral through the ABL: 
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radiation up from surface 

radiation into surface 

radiation up from ABL 
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latent heat of evaporation 
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Measurements or estimates from the CEAREX drift, during periods 
when reliable data were available, were used to derive ranges of the values for each 
of the terms in Equation (5) and are given in Table 7. The time derivative (Term 
1), advection (Term 2) and cloud top cooling (Terms 7 and 8) dominate the mean 
temperature equation. Surface fluxes (Term 3) and entrainment (Term 4) are an 
order of magnitude smaller and clear air radiation divergence (Terms 5, 6, 7 and 8) 
and latent heat (Terms 9 and 10) are almost another two orders of magnitude 
smaller. The smallness of the latent heat term justifies ignoring latent heat when 
considering the effect of clouds on ABL temperature. This is based on periods 
when ABL temperatures are below -20 C. As shown by Curry (1986), the latent 
heat term can be important in the summer when ABL temperatures are near 0 C. 

3. Graphical Representation of ABL Temperature and Surface Heat Flux 
Events. 

Only terms 1 and 3 in Equation (5) were directly measured; the others 
were estimated for typical Arctic conditions. A plot of the measured terms, Term 1 
vs. Term 3, provides a visualization of factors affecting ABL temperature 
(diagrammed in Figure 18). The vertical axis represents surface turbulent heat flux 
while the horizontal axis represents time change in ABL temperature. This 
represents a heat flux-temperature change or F-C diagram. A steady-state 
equilibrium (no turbulent surface heat flux) situation is represented by the origin in 
Figure 18. The dashed lines in Figure 18 represent lines of constant Term 2 + 
Term 4 through 9. These lines will be termed "isoadvects" because advection is 
usually the dominant term, although factors other than advection may be important 
for individual cases. 
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TABLE 7 

SCALE ANALYSIS OF ABL TEMPERATURE EQUATION 

(units are Wm’ 2 ) 



Term 


Percentile 

5 


95 


Median Magnitude 


1 . 


-121 


86 


34 


2. 


-110 


100 


35 


3. 


-12 


11 


4 


4. 


0 


10 


3 


5.+ 6. + 7. + 8. 
(Clear) 


-1 


1 


<1 


5. + 6. + 7. + 8. 
(ABL Cloud) 


-90 


0 


35 


9.+ 10. 


-3 


3 


1 
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Figure 18 F-C Diagram Showing Characteristics of Different Types of Sensible 
Heat Flux Events. The X-Axis represents the local time derivative of temperature 
times ABL depth and volumetric heat capacity (Wm' 2 ) and the Y-axis represents 
sensible heat flux (Wm’ 2 ). Pure Type R events fall along line AB while pure Type 
A events fall along line segments CD or EF. 
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Pure Type R situations (i.e. Terms 1 = Term 3) would have points only 
along the zero isoadvect, which is represented by line AB in Figure 18. A cloud- 
clearing event occurring in a situation which was previously in equilibrium would 
cause the representative location in Figure 18 to move from the origin to Point A as 
the negative surface fluxes respond to the radiation change. Eventually, the cooling 
of the ABL will counteract the fluxes caused by the new radiation conditions and the 
representative location will move back to the origin. Similarly, a cloud-forming 
event will move the representative location toward Point B and then back to the 
origin. Note that for pure Type R situations, the cloud-clearing or cloud forming 
must occur above the ABL; clouds within the ABL can affect ABL temperature 
directly, a Type A process. 

A pure Type A situation (surface temperature driven entirely by the 
turbulent heat flux) would begin at a point on the X-axis because at first the ABL 
temperature will not have changed and the fluxes will be zero. For example, at the 
onset of steady cold advection, the representative location would be on the X-axis at 
the value of the advection, such as Point C. Eventually the cooling ABL will 
generate positive fluxes and the isoadvect will be followed until Point D is reached. 
At this time, the ABL has cooled enough so that the positive heat flux is balanced by 
the advection and a steady-state is reached. Alternatively, a warm advection event 
would cause the representative location to move from Point E to Point F. Note that 
line AB is entirely within the upper right or lower left quadrants, where the signs of 
Term 1 and Term 3 are the same. Lines'CD and EF are in the upper left and lower 
right quadrants where the signs of Term 1 and Term 3 are opposite. This is 
consistent with the earlier, more general definitions of Type R and Type A 
situations which were defined in Table 6. "Pure" Type R and "pure" Type A events 
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were shown in Figure 18. In nature, completely "pure" situations almost never occur 
and advection is not constant. The more general definitions for Type A or Type R 
categories will be used for this study. Therefore any situation with changing 
temperatures and surface heat fluxes must be either Type A or Type R. 

F. EXAMPLES FROM CEAREX DRIFT 

Term 1 vs. Term 3 from Equation (5) for the entire CEAREX drift period 
when surface flux measurements were available is plotted in Figure 19. The data 
are averaged over one hour intervals. The vertical scale in Figure 19 is exaggerated 
compared to Figure 18 so that the isoadvects are almost vertical. The majority of 
the points are clustered toward the center where fluxes and temperature changes 
are small. There is too much noise (randomness) in the data to detect any 
discernable patterns in the central data cluster. 

The large heat flux events (high y-axis magnitude points in Figure 19) are 
examined more closely because these have the potential to strongly influence the 
effect of atmospheric stability on wind stress. Unlike the center of the scatterplot, 
the outer points in Figure 19 show some definite trends. The largest negative heat 
fluxes (below -15 Wm' 2 ) usually occur during periods of dropping temperatures. 
This is consistent with a cloud clearing event which cools the ABL by first cooling 
the surface. Type R events are more common than Type A events for explaining 
large negative heat fluxes. This study will concentrate on Type R events, since these 
are the situations when clouds have the strongest effects on wind stress. 

A good example of a Type R event occurred from 1800 8 October to 1200 9 
October (Figure 20). At the beginning of this period, the representative locations 
on the F-C diagram were near the origin, but then they moved down and to the left. 
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Sensible Heat Flux (Um A 2) 




Figure 19 F-C Diagram of Entire CEAREX Drift. Heat flux scale exaggerated 
compared to Figure 16. 
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Sensible Heat Flux (Um~2) 




Figure 20 F-C Diagram of Clearing Event From CEAREX Drift 1800 8 Oct to 1200 
9 October. The axes scaling is different compared to Figure 18. 
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Longwave radiation and sky observations (not shown) clearly indicate that the 
observed decrease in temperature and in negative sensible heat flux was associated 
with a change in cloud conditions. This is a Type R case since the points are near 
the zero isoadvect. 

Most points in the lower left quadrant of Figure 19 are to the left of the zero 
advect. This means that the ABL cooling cannot be entirely explained by a one- 
dimensional, Type R, surface heat flux argument for most cases. The one- 
dimensional reasoning would require that cloud clearing occur everywhere at once. 
This does not happen; the cloud edge moves horizontally, usually in the same 
direction as the air parcels in the ABL. Therefore, the cooling that occurs after 
cloud clearing is more intense because the air parcels have been exposed to the 
clear conditions and cold surface for a longer period than the time after the cloud 
clearing event at any particular location on the surface. 

Large upward (positive) sensible heat fluxes are equally divided between Type 
A and Type R events (top of Figure 19 to the left and right respectively). Cloud- 
forming events do not seem to dominate the positive upward heat flux cases to the 
extent that cloud clearing dominates negative heat flux events. 

The largest temperature increases (to the right of 150 Wm’ 2 in Figure 19) are 
all associated with negative fluxes, a Type A situation. But large temperature 
increases (to the left of -150 Wm’ 2 ) can occur with both upward and downward heat 
fluxes. Again, this indicates that cloud-clearing events are a relatively more 
important cause for surface heat fluxes than cloud-forming events. This means that 
the magnitude of the ABL warming associated with cloud formation is less than the 
magnitude of the ABL cooling associated with cloud clearing. 
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Two reasons explain the difference in surface heat flux response between 
cloud-forming and cloud-clearing situations. The first reason is that cloud-forming 
events often involve ABL clouds. The radiational warming is counteracted by cloud 
top cooling so that there are large positive fluxes at the surface but the ABL 
temperature does not warm. Cloud-clearing has no counteracting radiational 
warming in the ABL. The second reason is that the depth of the ABL is usually 
lower for cloud-clearing than cloud-forming events. Shallow ABLs respond more 
quickly to surface heat fluxes than do deep ABLs. 

Several conclusions concerning the effect of clouds on ABL processes can be 
derived from the CEAREX drift data discussed in this sub-section. Type R ABL 
events occurred at about the same rate as type A events when all data are 
considered. Therefore, cloud condition changes appear to be as important as other 
factors in affecting surface heat fluxes when all cases are considered. Large 
negative heat fluxes are particularly well-correlated with cloud condition changes. 

This section demonstrates that clouds are often important to surface heat 
fluxes. Therefore, much of the wind stress variation caused by surface heat fluxes 
over sea ice or ABL stability effects can be attributed to cloud condition changes. 

G. CONCLUSIONS ON CLOUD EFFECTS ON ATMOSPHERIC 

THERMODYNAMIC STRUCTURE 

This examination shows the strong association between cloud conditions and 

the thermodynamic structure of the lower atmosphere over sea ice in winter. 

1 

Observed data show that clouds have a large effect on surface temperature which in 
turn is closely linked to the ABL properties by turbulent fluxes. Clouds prevent the 
formation of surface-based inversions and are associated with higher inversion bases 
and multiple-level inversions. The presence of low clouds leads to no correlation 



66 



between the inversion height and wind speed. During the CEAREX drift, clouds 
had as great an effect on surface layer stability as all other factors combined. Cloud 
clearing events were especially effective at generating relatively large negative 
surface heat fluxes and cooling the ABL. 
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IV. THE EFFECT OF CLOUDS ON SURFACE LAYER STABILITY, THE 
SNOW/ICE LAYER AND WIND STRESS 



A. DEFINING SURFACE LAYER STABILITY EFFECTS 



The quadratic geostrophic drag coefficient, C G 2 , parameterizes the effect of 

'y 

surface and ABL physics on wind stress. The physical processes which affect C G * 



can be divided into two categories (1) surface layer effects and (2) outer ABL layer 
effects. The quadratic geostrophic drag coefficient can be expressed as 



where U 10 is the wind speed 10 meters above the surface, which is near the top of 

the surface layer. Surface layer effects on wind stress are parameterized by the 
surface drag coefficient, C d , while outer ABL processes are parameterized by the 
quadratic reduction factor, C UG . 

Surface heat fluxes cause surface layer stability effects which change the value 
of C d . The changes in C d caused by clouds will be analyzed in this chapter. The next 

chapter will include the effects of clouds on the reduction factor. 

B. PURPOSE OF SURFACE LAYER STABILITY STUDIES 

Discussion in Chapter III addressed how changes in cloud conditions caused 
heat fluxes and stratification at the surface. Here, these heat fluxes will be 




( 6 ) 
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analytically and numerically modeled assuming pure Type R events, i.e. no 
advection, entrainment or diabatic ABL cooling other than surface fluxes. One 
purpose of this section is to determine analytical expressions for important time 
scales related to Type R heat flux events. Another purpose is to determine the most 
important parameters for predicting heat flux and wind stress on various time scales 
over sea ice. 

OG modeled the effect of an instantaneous change in downward radiation on 
T air and surface heat flux (Figure 21) to describe changes in conditions above the 

ABL from cloudy to clear and back to cloudy. They did not consider the wind stress 
implications of the radiational changes. Heat flux values from Figure 21 indicated 
that there were periods (the first two hours of day 2) when turbulence and wind 
stress would be almost completely suppressed and other periods (first hour of day 6) 
when it would be enhanced by 30%. These conclusions were based on surface layer 
Monin-Obukhov similarity theory (MO). Inclusion of ABL changes would magnify 
the surface layer stability effects. Although the OG case (Figure 21) had low wind 
speeds and therefore may not be significant to wind stress variations, it shows that 
for certain conditions, changes in cloud cover can have large relative effects on wind 
stress for short time periods. 

C. CONDUCTION OF HEAT THROUGH SNOW/ICE (ANALYTICAL 

MODEL) 

A skin depth, d skin , is defined as the effective or e-folding distance that 

longwave radiation penetrates the snow/ice surface. If there is an instantaneous 
change in downward surface radiation, F rad , from a previously equilibrium (no net 
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TIHEtORYS) 



Figure 21 Model Simulation Time Series of (a) Surface Temperature, (b) 
Difference between Surface and 10 meter Air Temperature, (c) Surface Sensible 
Heat Flux. A change from cloudy to clear occurs at 2 days and becomes overcast 
again after 6 days. From Overland and Guest (1991). 
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flux) situation, then the surface temperature, T sfc , cannot change temperature any 
faster than 



djsic 

dX 



rad 



^skin^skin ’ 



( 7 ) 



where C skin is the volumetric heat capacity of the snow or ice at the surface. 

The heat conduction within the snow or ice cover is governed by the one- 
dimensional heat flux equation (Semtner, 1976) 

3T dh 

C— = kT , (8) 

<9t dz 



where k is the thermal conductivity of the snow or ice and C is the volumetric heat 
capacity. In the first few seconds following a change in radiation conditions, the 
surface temperature has not changed enough to generate sensible turbulent fluxes. 
Therefore, all the radiational imbalance not stored in the skin layer must be 
counteracted by a conductive flux, F c , from the skin layer to the snow just below. 

F„<, = F c = k*(f ) sfc (9) 



If it is assumed that C and k are constant to infinite depth and d s ^j n is 
infinitesimal, the only parameters governing the value of T sfc are F rad (Js’W 2 ), C 
(Jm^K* 1 ), k (Jm'V'K* 1 ) and the time since the change in radiation conditions, t (s). 
From dimensional analysis 
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( 10 ) 



T S fc(t) = T sfc (0) - nF rad 




The dimensionless constant (n) was determined to be equal to 1.13 based on 
results from a numerical model similar to Semtner (1976), with 1000 computational 
levels in the snow layer and a time step of 0.01 seconds. 

A time scale, t skin , can be derived from Equations (7) and (10) which 

represents the time required to heat the skin layer enough so that all the radiational 

5T sfc 

heat imbalance is conducted to the lower layers and ^ no longer depends on 
storage of heat in the skin layer: 



t 



skin 



(nd skln ) 2 

^skin 



( 11 ) 



where v snow = ^ is the thermal diffusivity. At approximately t < t sldn Equation (7) 

is applicable while at t > t skin Equation (10) applies. Note that t skin is independent 
of F rad . Using a value of 1 mm for d skin and v = 1.0 x 10 which are characteristic 
values for feathery new snow (Stull, 1988), t skin has the value of 12.8 seconds. For an 
ice surface, t skjn equals 1.4 seconds. Old snow has intermediate values. In nature, 
downward longwave radiation changes usually takes a few minutes; therefore the 
heat storage represented by Equation (7) is not significant for downward radiation. 
Equation (7) would be more applicable to downward solar shortwave radiation, 
which is not present in the winter. Shortwave radiation can change within a few 
seconds as the sun goes behind a cloud and penetrates several centimeters into the 
snow/ice. 



72 



Equation (10) is not valid after the time, t turb , when T sfc has changed enough so 
that the surface turbulent heat flux, F^, is significant, assumed to be when the 



sensible heat flux is 10% of the initial radiational imbalance 



F, 



sens 




( 12 ) 



The turbulent heat flux can be approximated by 



F 



sens 



ATp CpU 10 C] 



10 '-H » 



(13) 



where p is the density of air, c p is the specific heat of air at constant pressure, aT is 
the difference in potential temperature between the surface, T sfc , and a reference 
height near the top of the surface layer, T 10 . U 10 is the wind speed at the reference 
height, and C H is the heat transfer coefficient. Assuming that Equation (10) is valid 
and T 10 remains equal to T sfc (t = 0) when t < = t turb , an approximate expression 
for t turb can be derived using Equations (10), (12) and (13): 



This time scale represents the time required for the turbulent heat fluxes to become 
significant after a change in radiation conditions. Some values for various snow/ice 
types are listed in Table 8. The low valties of t turb explain the previously described 

measurements of the close coupling between the turbulent heat fluxes and radiation. 



Ck 



t 



turb 



(10 p c p Ujo C H n) 



(14) 
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TABLE 8 



VALUES OF T turb (seconds) FOR U 10 = 5 ms' 1 , C H = 1.0 x 10' 3 



New feathery snow (Stull, 1988) 


0.7 


Old packed snow (Stull, 1988) 


44.8 


Ice (Stull, 1988) 


614.0 


Snow (Semtner, 1976) 


33.9 


Ice (Semtner, 1976) 


604.0 


CEAREX drift top 3 cm (estimated) 


20.0 
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E. CONDUCTION OF HEAT THROUGH SNOW/ICE (NUMERICAL 

MODEL) 

1. Model Description 

The analytical models described above are useful for determining simple 
formulae for time scales associated with radiation changes and surface temperature 
and can be accurate for predicting T sfc for the first few seconds or minutes. 
However, a numerical multi-level model is better suited for predicting T sfc when the 
radiational change is "felt" at lower levels in the ice/snow, where k and C may not 
be constant. 

The numerical snow/ice thermodynamic model used was based on 
Semtner (1976). This model numerically solves Equation (8) for each level. The 
surface temperature, T sfc , is determined by assuming a balance of net radiation, 

turbulent sensible and latent heat fluxes, and heat conduction through the ice at the 
surface. The value of T sfc for the next time step is based on analytical expressions 
for the partial derivatives of each of the fluxes with respect to a change in T sf , 
Semtner (1976) prescribed the turbulent fluxes since he was studying long term 
processes associated with ice growth and climate. For this study, time dependent 
expressions for the effect of turbulent fluxes on T S f C have been added. 

2. Model Time Scales 

The purpose of this sub-section is to provide a guide for choosing enough 
vertical layers so that turbulent fluxes can be accurately modeled. Appropriate time 
scales associated with various snow/ice model parameters are described. These 
time scales can be used to determine how many layers are needed in the snow/ice 
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